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Abstract: We analyze the effect of phase fluctuations in an optical 

communication scheme based on collective detection of sequences of 
binary coherent state symbols using linear optics and photon counting. 
When the phase noise is absent, the scheme offers qualitatively improved 
nonlinear scaling of the spectral efficiency with the mean photon number 
in the low-power regime compared to individual detection. We show that 
this feature, providing a demonstration of superaddivitity of accessible in¬ 
formation in classical communication over quantum channels, is preserved 
if random phases imprinted on transmitted symbols fluctuate around a 
reference fixed over the sequence length. 
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1. Introduction 

Binary phase shift keying (BPSK) is a well established modulation format in optical comu- 
nication IIIlEl. It employs a pair of coherent states with the same mean photon number n and 
opposite phases + = e'® and — =e‘^. In the low-power regime, when n<^l, its readout becomes 
non-trivial because of quantum fluctuations, which prevent perfect discrimination between the 
two symbols constituting the BPSK alphabet. For individual detection of BPSK symbols, the 
minimum discrimination error is given by the celebrated Helstrom bound and can be attained 
using the Dolinar receiver Il3l4l based on photon counting and fast feed-forward operations. The 
resulting spectral efficiency, given by Shannon mutual information expressed in bits, scales in 
the low-power regime linearly with the mean photon number as /hbI ~ (2/ln2)H. Quantum 
mechanics offers a possibility to change qualitatively the above scaling through collective de¬ 
tection of multiple symbols. The ultimate limit of such a general strategy is given by the Holevo 
quantity x 0 which in the case of BPSK expands for n<^l as x ^ exhibiting en¬ 

hanced nonlinear scaling with the mean photon number. The ability to improve the attainable 
transmission rate by implementing a joint measurement on sequences of symbols is referred to 
in quantum communication theory as the superadditivity of accessible information 0. 

Recently, a linear-optics scheme to attain superadditivity for the BPSK format has been pro¬ 
posed by Guha 0. The basic idea is to prepare sequences of BPSK symbols defined by columns 
of a Hadamard matrix. Because of orthogonality properties of Hadamard matrices, applying an 
appropriate linear-optics transformation to such Hadamard words maps them onto the pulse po¬ 
sition modulation (PPM) format that is suitable for direct detection. This strategy can approach 
in the leading order the Holevo quantity when optimized over the sequence length 0|9l. The 
purpose of this paper is to study the impact of excess phase noise on collective BPSK with 
Hadamard words. As the main result, we demonstrate that the nonlinear scaling of spectral ef¬ 
ficiency in the low-power regime is retained when individual BPSK symbols experience small 
phase fluctuations occurring around a reference locked over the sequence length. Remarkably, 
the enhanced scaling persists despite the optimal length of Hadamard words growing indefi¬ 
nitely large as the mean photon number tends to zero. 

This paper is organized as follows. In Sec.|2]we review collective BPSK with Hadamard 
words and introduce the phase noise model. Optimization of the spectral efficiency is consid¬ 
ered in Sec.0 Finally, Sec.|4]concludes the paper. 

2. Collective BPSK with Hadamard words 

Hadamard matrices are square symmetric matrices with entries hp = hij = ±1 and mutually 
orthogonal columns, or equivalently rows. They exist for dimensions L equal to powers of 
2 and selected other integers. In the proposed collective BPSK scheme 0 shown in Fig. [T] 
sequences of input symbols, called in the following Hadamard words, are chosen as columns 
of a Hadamard matrix, i.e. the /th word is given by {h\i,h 2 iT ■ ■ ,hu), where I = 

These symbols define the complex amplitudes of individual pulses, i.e. the amplitude of the jth 
pulse in the /th word is hjis/^, where h is the mean photon number per time bin. 

At the detection stage the sequence, after synchronizing the time bins, is fed into a linear 
optics circuit that implements a unitary transformation of input complex amplitudes described 
by a rescaled Hadamard matrix {h^j/s/L). The output amplitude at the kth port is thus given 
by \/n/LY,jhkjhji = \/Tli5ki, where in the second expression we have made use of the or¬ 
thogonality properties of Hadamard matrices and denotes Kronecker delta. It is seen that 




Fig. 1. Sequences of BPSK symbols ± shown on the left are prepared as rows (or equiva¬ 
lently columns) of a symmetric Fladamard matrix. A linear circuit, described by a rescaled 
Fladamard matrix, transforms them into the PPM format visualized as tall solid red pulses 
localized in single output bins. When sequences are subject to phase noise, some of the 
input sequence energy becomes evenly distributed across all other bins, depicted with light 
red pulses. Assuming no dark counts, the sequence is either identified correctly with prob¬ 
ability p by the position of the detector click (solid black arrows) or erased when no pho¬ 
tocounts are generated (dashed blue arrows). In the presence of phase noise, the fraction of 
the sequence energy spread over all remaining PPM bins may produce a click in a given 
wrong bin with a probability q (dotted purple arrows). Events when clicks occur in two or 
more bins are treated as erasures. 


amplitudes are zeroed at all output ports of the circuit except the (th one, which carries the 
entire sequence energy characterized by the mean photon number Ln. The non-empty port, and 
consequently the input Hadamard sequence, can be identified using photodetectors monitoring 
individual output ports. Neglecting dark counts and assuming that binary on/off detectors are 
used, identification is unambiguous when one of the detectors generates a photocount. In the 
absence of any detection event, no information about the input sequence is recovered. 

From the communication theory viewpoint the above scheme is described by an erasure 
channel nni with L equiprobable values of the input variable that is either transmitted faithfully 
through the channel or replaced by an erasure flag. In our case, the erasure probability is given 
by exp(—L«), which corresponds to registering no photocount at any output port. For such 
a scenario. Shannon mutual information per time bin reads T^^[l — exp(—L«)] log 2 L. When 
n 1, the optimal sequence length can be found approximately by expanding 1 — exp(—L«) Ri 
Lh — [Lnf'/2 mm. This calculation yields spectral efficiency with the same leading-order 
dependence on h as the Holevo quantity X log 2 (l /n). 

Suppose now that the transmission of each BPSK symbol is affected by uncorrelated Gaus¬ 
sian phase noise im characterized by a standard deviation a. In this model, before the readout 
stage the phases of individual BPSK symbols are (e'‘^*/ 2 i/,e'‘^/ 22 /, ■ ■ ■ where L random 

phase variables , 02,..., are described by the same Gaussian distribution 

p(0) = exp[—0^/(2(7^)]/V27r(J^. (1) 

Because of the phase noise, at the output of the Hadamard circuit the energy of the sequence 
will no longer be concentrated at a single output port. Consequently, the communication scheme 
cannot be represented as a standard erasure channel, but a possibility of incorrect identification 
of the input sequence also needs to be taken into account. 

Let us analyze the above more general scenario in quantitative terms. For the /th Hadamard 
word and an individual realization of all random phase variables, the kth output port of the 


















linear optics circuit carries a coherent state with an average photon number 
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We will be interested only in events when exactly one detector clicks. Events when two or 
more detectors produce photocounts will be treated as erasures. The exclusive probability for 
the detector monitoring the kth port to generate at least one photocount, while all other do not 
click at all, reads 


" l)> (3) 

where the angular brackets denote the statistical average over the phase noise. The last sim¬ 
plified expression for in Eq. Q is obtained from energy conservation in the course of the 
linear circuit transformation, which implies that Lj=i = Lfi. It is easy to see that the proba¬ 
bility assumes one of only two values depending on whether k = l or k^l. This is because 
in Eq. (|2|i all the products hi^jhji are equal to one when k = I and they are split evenly between 
1 and — 1 when k^l, the latter property following from the orthogonality of Hadamard matrix 
columns. Consequently, it is sufficient to introduce two parameters; the probability of correctly 
identifying the input Hadamard word p = p^p and the probability of obtaining a wrong answer, 

q = pP for k ^ 1. The resulting communication scheme can be viewed as a generalized era¬ 
sure channel with uniform noise also depicted in Eig. [T] which in addition to deleting the input 
variable can also scramble its value. 

The formulas for the parameters p and q following from Eq. Q involve arbitrarily high mo¬ 
ments of the distributions of random phase variables and are difficult to handle analytically. 
One can however derive simple bounds on p and q that will prove useful when optimizing the 
sequence length. The starting point is to average over the phase noise the mean photon number 
at individual output ports. This yields 

{liP)=n[\ + {L5u-\)e-^\ (4) 

(') ; (')\ 

In the case of p, the convexity of the exponential function implies that {e^i ) > K Using 
Eq. dUi for k = / yields 

p > exp[—(L—l)n(l —)] —exp(—Ln). (5) 


On the other hand, the exclusive probability of a click is always less or equal than the inclusive 
one, which for k ^ I implies that 

q<{\-e^^kp<(^pP)=n{\-e^’^). ( 6 ) 

The above bounds are conservative, underestimating the probability of a correct identification 
of the Hadamard sequence while overestimating the probability of an identification error. 


3. Optimization 


Starting from the standard definition ifTOll . a lengthy but straightforward calculation shows that 
Shannon mutual information per time bin for a generalized erasure channel with uniform noise 
can be written as 


p 1 , 

I = j[p- 
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Fig. 2. (a) Shannon mutual information / given in Eq. 0 as a function of the sequence 
length L for h = 10^^ and several dephasings CT calculated for p and q estimated in Eqs. (0 
and (solid lines) and using Eq. 0 expanded up to (dashed lines). The case cr = 1.7 
illustrates the breakdown of approximations used in the analysis, (b) The ratio /(a)//o 
characterizing the effect of phase fluctuations on mutual information as a function of the 
average photon number h. Information /(cr) is optimized over the sequence length L treated 
as a continuous parameter taking bounds given in Eq. 0 and 0 (solid lines), numerical 
expansions up to (dashed lines), and the closed approximate formula from Eq. 0 (dotted 
lines). The reference value Iq has been obtained by optimizing numerically over L the exact 
expression for mutual information in the noiseless case. Arrows indicate asymptotic values 
For CT = 0.3 dashed and solid lines overlap. 


where H{x) = —.xlog 2 Ji: — (1 — x) log 2 (l — x) denotes binary entropy. In Fig. |2]we present I as 
a function of the sequence length L for n = 10^^ and increasing dephasing calculated using 
two approaches. The first approach is to compute numerically the probabilities using Eq. (0 
expanded up to the second order in n, while the second one uses bounds derived in Eqs. 0) and 
0. It is seen that both calculations yield very close locations of the optimal sequence length L. 
Eurthermore, the optimal L does not change noticeably with the strength of dephasing as long as 
the latter does not become too large. One should also note that the permitted values of L defined 
by the dimensionality of Hadamard matrices are evenly distributed in the logarithmic scale 
used for the abscissa in Eig. |2ja). Eollowing this observation, we will treat L as a continuous 
parameter when optimizing mutual information. 

Let us now inspect closer the expression for mutual information given in Eq. (|7]i with the 
aim to derive an approximate analytical formula for its optimal value in the regime of weak 
dephasing. Motivated by the numerical example shown in Eig. |2a), we will assume that L ^ 1. 
The self-consistency of this assumption will be confirmed afterwards. It is easy to verify that 
the factor [p -F (L— \ )q]/L multiplying the binary entropy in the second term is equal to h in 
the leading order, whereas the last term in Eq. (|2ll is of the order OinjV) and therefore can be 
neglected for long sequences. Eurther, applying a linear expansion in h to p and q in the second 
term of Eq. (|7]i, using Eq. (|4]i, and taking L ^ 1, mutual information can be simplified to an 
approximate form I r; {p/L) log 2 L —nH{e^^ ). Explicit dependence on L occurs only in the 
first term, which has a form analogous to mutual information for a standard erasure channel 
with the non-erasure probability equal to p. 

Given this analogy, we will use the approach presented in |0| and expand p up to the second 
order in Ln: 

p^e-^'Ln-^j{2e^^-\){e-^^Lnf. ( 8 ) 

This makes our problem equivalent to optimizing the PPM format when the average photon 
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number per bin is equal to n and the non-zero pulses exhibit super-Poissonian photon 











statistics characterized by the normalized second-order intensity correlation function = 

2 __ 

26°^* — 1. The solution to this problem found in ||9J satisfies assumptions made about the optimal 
sequence length Lopt, i.e. n Topt« ^ 1- The approximate expression for mutual information 
based on this result can be compactly written as 

1(a) ft! e-‘^"Hn((2 - e-^^)n) - nH(e-^'), (9) 

where n(x) is the photon information efficiency function for the PPM format with Poissonian 
pulse photon statistics given explicitly in ||9l- 

Figure |2b) depicts the ratio 1(a)/I q of mutual information optimized over the sequence 
length in the presence of phase noise 1(a) to the ideal noiseless case /o as a function of 
the average photon number n. Two approaches used to produce Fig. |3a) are found in good 
agreement with the closed analytical expression derived in Eq. (|9]l, which confirms the con¬ 
sistency of our analysis. Using in Eq. (|9]l the asymptotic form of the photon information ef- 
hciency function n(x) ~ log 2 (l/x) for small arguments x <C 1 yields in the leading order 
1(a) ~ hlog 2 (l/h). Consequently, the ratio I(a)/Io tends for « —0 to a constant value 
, marked with arrows in the ordinate of Eig. |2jb). It is seen that this asymptotic value is 
attained faster with decreasing h for weaker dephasings a. Importantly, the Holevo-type scal¬ 
ing of mutual information in the presence of phase noise is retained and the principal effect 
of phase fluctuations in the asymptotic regime fi ^ 0 is introduction of a multiplicative factor 

4. Conclusions 

The communication scheme based on Hadamard words constructed from BPSK symbols with 
a linear-optics joint-detection receiver Q attains superadditive accessible information through 
collective measurements of transmitted optical signals. Quantum-enhanced detection schemes 
are potentially more sensitive to noise and imperfections, as recently noted in the case of indi¬ 
vidual quantum receivers for the BPSK format ifTTIl . Theoretical analysis presented in this work 
shows that a moderate amount of phase noise can be tolerated in the collective BPSK scheme 
without compromising the Holevo-type scaling of mutual information with the mean photon 
number. 

In our calculations, the phase variables fluctuated around a reference value that remained 
constant over the entire Hadamard sequence. An alternative model to analyze would be free- 
floating Markovian phase diffusion discussed recently in ifT^ for a different keying scheme. It 
should also be noted that there exist coherent strategies capable of reading out the PPM format 
more efficiently compared to direct detection Ha. More generally, it would be interesting to 
investigate practical schemes to attain collective enhancement in the case of larger alphabets 
such as quadrature phase shift keying for which individual quantum receivers operating below 
the shot noise limit have been recently demonstrated inaiia. 
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